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Abstract 

The concentrations of 16 US-EPA criteria polycyclic aromatic hydrocarbons (PAHs) were monitored during various 
seasons at six different cities/locations in the northern part of Belgium. Pressurized liquid extraction (PLE) was used for 
the fast recovery of PAHs from quartz fibre filters (QFFs) and polyurethane foams (PUFs) in less than 30 min with 
minimum solvent consumption prior to their analysis with high-performance liquid chromatography. The concentrations 
of PAHs varied significantly at the studied sites and showed relation to different anthropogenic activities, such as vehicular 
emission (diesel/gasoline), incinerator, petroleum/oil burning, coke production, and wood/coal combustion. The annual 

_o 

average PAH concentration ranged from 17 (at a rural site) to 114ngm (near a petroleum harbour and industry). Most 
of the human carcinogenic PAHs were found to be associated with suspended particulate matter, which forms around 
~55% of the total PAH levels in aerosols. Relatively higher concentrations of PAHs were detected in aerosol samples 
during winter compared with other seasons, whereas no clear seasonal trend was observed for the vapour phase PAHs. 
This fraction is likely to be more local in origin; hence, it can be used as a site-specific characteristic. The site-specific 
concentrations of individual PAHs were also used in diagnostic ratio evaluations and in principal component analysis to 
find their relation with various anthropogenic activities. These results show that the vehicular emission is a major source of 
PAHs in Flanders, although other anthropogenic sources, as above, have also an impact on the total PAH levels. 

© 2005 Elsevier Ftd. All rights reserved. 
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1. Introduction 

Polycyclic aromatic hydrocarbons (PAHs) are of 
major health concern, mainly due to their well- 
known carcinogenic and mutagenic properties 
(Ravindra et al., 2001). PAHs and their metabolites 
can form DNA adducts, which can induce muta¬ 
tion. These adverse properties, together with the still 
growing presence of PAHs in the environment, 
make their monitoring to be a top priority task in all 
programs for the evaluation of both the human 
health risk and the environmental hazards. 

As it is known, PAHs are products of incomplete 
combustion and pyrolysis of fossil fuels and other 
organic materials from natural and anthropogenic 
sources. In urban and industrial atmospheres, 
PAHs are almost entirely of anthropogenic origin 
(Jones et al., 1989). PAHs are present in the 
atmosphere both in the vapour and the condensed 
(aerosol) phases (Bidleman et al., 1986). Low 
molecular weight PAHs tend to be more concen¬ 
trated in the vapour phase, while the ones with 
higher molecular weight are often associated with 
particulates. Thus, it is important to reveal the 
abundance, distribution and potential sources of 
PAHs in the ambient air both in aerosols and in the 
vapour phase, so that air pollution caused by PAHs 
is to be efficiently controlled. 

The main sources of PAHs in towns and cities are 
road traffic and certain industrial processes (Lee et 
al., 1995; Ravindra et al., 2001). In order to enhance 
the accuracy of emission source identification, 
principal component analysis (PCA) was used in 
many studies (Ho et al., 2002; Park et al., 2002; 
Fang et al., 2004; Dallarosa et al., 2005). Addition¬ 
ally, as the PAHs associated with suspended 
particulate matter (SPM) vary significantly as a 
function of their emission sources, some PAH 


concentration ratios have been used to indicate 
vehicular emission sources (Rogge et al., 1993; 
Venkataraman et al., 1994; Khalili et al., 1995). The 
amount and type of PAHs in urban and rural 
areas exhibit seasonal fluctuations with generally 
higher concentrations in the winter season (Cotham 
and Bidleman, 1995; Harrison et al., 1996; Ravindra 
et al., 2006). 

In the present study the concentration of PAHs 
has been evaluated in Flanders, the northern part of 
Belgium, with the application of a fast analytical 
approach for their extraction. Samples were col¬ 
lected from the ambient air in both the vapour and 
particulate phases and the variation in PAH 
concentrations were assessed during different me¬ 
teorological and seasonal conditions. A further 
purpose of the study was to identify the site-specific 
characteristics with relation to PAH concentrations 
and to apply PCA to find the possible emission 
sources of these compounds. 

2. Experimental 

2.1. Sampling sites 

Six sampling sites were selected on the basis of 
different anthropogenic activities in Flanders, 
Belgium (Fig. 1). These sites include the following 
locations: Wingene, Hasselt, Mechelen, Antwerp- 
Borgerhout, Antwerp-Petroleumkaai and Zelzate. 
Two sampling campaigns were performed at each 
site during different seasons. From each campaign, 
7-9 samples were stochastically collected. Table 1 
shows the sampling period of each sampling site 
and the season(s) during that period. The most 
important meteorological parameters, such as the 
wind speed (WS) and wind direction (WD), the 
relative humidity (RH), the air temperature (7), 



Fig. 1. Locations of various sampling sites in Flanders (Belgium). 
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Table 1 

Sampling periods during the monitoring of PAHs 


Campaign/location 

Sampling periods 

Campaign 1 

Petroleumkaai (A) 
Borgerhout (A-W) 
Zelzate (W) 

Hasselt (W-Sp) 
Wingene (Sp) 

Mechelen (Sp-Su) 

22/09/01-09/10/01 

24/11/01-06/12/01 

20/12/01-18/01/02 

12/02/02-22/03/02 

06/04/02-14/05/02 

24/05/02-26/06/02 

Campaign 2 

Wingene (Su) 

Zelzate (Su-A) 

Hasselt (A) 

Mechelen (A-W) 
Petroleumkaai (W) 
Borgerhout (W-Sp) 

27/06/02-06/08/02 

18/08/02-11/09/02 

06/10/02-04/11/02 

16/11/02-15/12/02 

04/01/03-02/02/03 

16/02/03-13/03/03 

(W—winter, Sp- 

-spring, Su—summer, A— 

-autumn). 


air pressure (AP) and the precipitation rate (PPT) 
were also monitored and recorded during the 
sampling campaigns near the corresponding sam¬ 
pling sites. 

2.2. Sampling site characteristics 

Sites selected for PAH measurements in Flanders 
were chosen on the basis of the local differences in 
the type, distribution and proximity of emission 
sources. Two of the sampling sites are located in 
urban areas, i.e., Borgerhout and Mechelen, while 
one of the sampling sites (Wingene) is situated in a 
rural area. Two other sites (Zelzate and Hasselt) are 
located in suburban areas. One site (Petroleumkaai) 
is of industrial nature, as it is situated near the 
petroleum harbour and several refineries of An¬ 
twerp. The Zelzate site is situated near a highway. 
There are considerable industrial activities nearby 
the sampling sites at Petroleumkaai (oil refineries 
and petrochemical plants), Zelzate (steel industry) 
and Mechelen (dye/paint). Samples were also 
collected near an area of high traffic intensity 
(Borgerhout), i.e., 40,000-50,000 cars per day. 

2.3. Sampling procedure 

For the determination of PAHs, atmospheric air 
samples were taken at 1.5 m height using polyur¬ 
ethane foam (PUF) high-volume samplers (Model 
GPS-1, Andersen Instruments, Village of Cleves, 
OH, USA). PAHs present in the particulate and 
vapour phases of the ambient air were trapped on 
quartz fibre filters (QFFs) (Gelman, Ann Arbor, 


MI, USA) and PUF plugs, respectively. The 
instrument was calibrated using a manometer and 
top-loading orifice plate (Graseby Anderson, OH, 
USA). An airflow rate of approximately 1801 min -1 
was used with a daily sampling time of around 
24 h. After sampling, each filter and PUF was 
separately wrapped in aluminium foils and stored in 
a freezer at —20 °C, until they were processed to 
sample preparation and analysis. The sampling 
schedule during the monitoring period is shown in 
Table 1. 

SPM and PM 2 5 fraction of ambient air were 
determined by the use of a four-stage Dekati 
impactor (Tampere, Finland). This unit was oper¬ 
ated at a height of around 1.7 m from the ground for 
the same 24-h period as the PUF sampler and with a 
flow rate of 301 min -1 . Nuclepore filters with a 
diameter of 25 mm and a pore size of 0.2 pm were 
placed onto each stage of the impactor. Moreover, a 
quartz filter of 47 mm diameter was also placed into 
the last “back-up” stage of the impactor to collect 
rest of the aerosols, i.e., the “back-up” fraction. The 
filters were weighed with a microbalance, before 
and after sampling, and the particulate fractions 
were calculated according to the weight difference. 
The SPM and PM 2 .5 were also calculated from 
the weights of the diverse fractions collected on the 
impactor stages. 

2.4. Sample preparation by pressurized liquid 
extraction (PLE) 

The samples were extracted using Accelerated 
Solvent Extractor (ASE, Dionex, Titan Way, CA, 
USA). As a first step of this procedure, the filters 
and PUFs were loaded into 34 and 100 ml extrac¬ 
tion cells, respectively. The extractions were per¬ 
formed using a pressure of 10 MPa at an oven 
temperature of 100 °C. A solvent mixture of 1:1 
dichloromethane/acetone was used for extraction. 
The flushing volume was 60% of the volume of the 
extraction cell. The full cycle was repeated three 
times, and after extraction, the instrument was 
purged at 1 MPa for 60 s with nitrogen. Pre-rinsed 
250 ml screw cap I Chem Vials were used to collect 
the extracts. Finally, before analysis, the extracts 
were enriched in Turbo Vap 500 (Zymark, Hopkin- 
ton, MA, USA) concentration workstation after 
adding 0.6 ml acetic acid to avoid evaporation till 
dryness and were redissolved in 1 ml acetonitrile. 
Blanks of the filters and PUFs were also extracted 
using the same procedure as aforementioned. 
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2.5. Analysis of extracted samples by HPLC 

The extracts were analysed by a high-perfor¬ 
mance liquid chromatography (HPLC) system 
(Waters, Milford, MA, USA) with programmed 
fluorescence detection (Perkin-Elmer LC240, Nor¬ 
walk, CT, USA) and UV detection (Waters) for 
acenaphthylene. Separation of the PAHs was 
accomplished using a Vydac 201TP (250 x 4.6 mm) 
column, with a gradient elution ranging from a 
50 + 50% acetonitrile-water mixture to 100% acet¬ 
onitrile in 20 min. The fluorescence of PAHs was 
monitored with an automatic adjustment of the 
wavelengths for each compound according to its 
retention time. The excitation wavelengths ranged 
between 260 and 300 nm, while the emission 
wavelengths laid between 380 and 465 nm. The 
wavelength of the UV detector was 325 nm. Linear 
calibration graphs were obtained for all the studied 
PAHs in the concentration ranges of up to 
80ngml _1 . PAH concentrations of less than a 
nanogram in the air samples could be quantified. 
The detected and quantified PAHs were: naphtha¬ 
lene (Nap), acenaphthylene (AcPy), acenaphthene 
(Acp), fluorene (Flu), anthracene (Ant), phenan- 
threne (PA), fluoranthene (FluT), pyrene (Pyr), 
benzo(a)anthracene (BaA), chrysene (Chr), ben- 
zo(b)fluoranthene (BbF), benzo(k)fluoranthene 
(BkF), benzo(a)pyrene (BaP), dibenzo(a,^anthra¬ 
cene (DahA), benzo(g,h,i)perylene (BghiP) and 
indeno(l,2,3-c,d)pyrene (Ind). For these PAHs, the 
limits of quantification (ngnU ), calculated for the 
volume of air sampled, were 0.034, 0.040, 0.047, 
0.020, 0.003, 0.007, 0.021, 0.004, 0.021, 0.016, 0.015, 
0.020, 0.017, 0.003, 0.003 and 0.010, respectively. 
The calculation of quantification limits and other 
details can be found elsewhere (Tirez, 1995). 

2.6. Quality assurance 

All analytical data were subjected to strict quality 
control. Blank and spiked samples were analysed 
with each set of samples. Spiked filter and PUF 
samples were determined with good precisions and 
recoveries. Recovery efficiencies were found be¬ 
tween 80% and 120% by using a certified reference 
material of the National Institute of Standards and 
Technology (NIST): SRM 1650a (Diesel Particulate 
Matter). The NIST standard PAH mixture (SRM 
1647d, Schmidt, Amsterdam, The Netherlands) is 
used for calibration of the analytical method and 
spiking of the test samples. The standardized test 


method for PAHs was accredited by Behest (http:// 
belac.fgov.be) under the number 216-T (norm EN- 
45001 ISO-9002). 

3. Results and discussion 

3.1. Fast extraction of PAHs with ASE 

The traditional extraction methods for solid 
matrices are based on Soxhlet- and ultrasonic 
extraction. Both of these methods are time- and 
labour-consuming and require large amounts of 
organic solvents. Therefore, new extraction methods 
have been developed, e.g., microwave assisted 
extraction, supercritical fluid extraction and appli¬ 
cation of PLE (Dean, 1998; Bowadt and 
Hawthorne, 1995; EPA, 1995; Richter et al., 
1996). Several studies have shown that these 
methods can be equally or even more efficient than 
Soxhlet extraction (David and Seiber, 1996; 
Brumley et al., 1998; Berset et al., 1999). 

Application of PLE has been shown for the 
recovery of PAHs from various environmental 
matrices such as contaminated soils (Lundstedt 
et al., 2000), solid waste (Popp et al., 1997), water 
samples (Doong et al., 2000), fly ash, and lignite 
coal (Arditsoglou et al., 2003) and for mussels 
(Martinez et al., 2004). Further, Godoi et al. (2004) 
have also shown the application of PLE for the 
extraction of PAHs from particulate matter col¬ 
lected on glass fibre filters followed by GC-MS 
analysis. However, to the best of our knowledge, 
there is not any study in the literature, which 
provides fast conditions for the recovery of PAHs 
from QFF and PUF samples collected from 
ambient air prior to HPLC analysis. In the 
experimental section, the optimum conditions for 
the maximum recovery of PAHs from QFF and 
PUF samples are provided. The results show that 
application of PLE offers very fast recovery of 
PAHs, i.e., in 28 min, with minimum consumption 
of toxic solvents. The total amounts of the solvent 
used for the extraction of QFF and PUF samples 
were 80 and 60 ml, respectively. These solvent 
volumes can be reduced further, either by the use 
of an extraction cell of appropriate size or by the 
use of inert material such as Hydromatrix (cleaned 
and sieved diatomaceous earth) as a solid support 
to fill the empty space in extraction cell, since 
the amount of solvent consumption depends on the 
sample size. 


K. Ravindra et al / Atmospheric Environment 40 (2006) 771-785 


775 


3.2. Evaluation of site-specific PAH concentrations 

The background levels of PAHs in the ambient 
air are reported to be 0.02-1.2 ngm~ in rural areas 
and 0.15-19.3 ngm~ in urban areas (Ravindra 
et al., 2001). Because of their low vapour pressure, 
PAHs in air are present at ambient temperatures 
either as gases or condensed-phase (liquid and/or 
solid) particles (i.e., SPM), or gases adsorbed on the 
particulate fraction. The site-specific concentrations 
of PAHs both in the vapour- and particulate 
fractions were determined at various sites in 
Flanders during different seasons. As can be seen 
in Table 2, the highest annual average concentra¬ 
tion of PAHs was found at Petroleumkaai 
(114ngm - ), followed with lower concentrations 
by Borgerhout (54ngm _ ), Zelzate (52ngm _ ), 
Hasselt (43ngm - ), Mechelen (42ngm~ ) and 
Wingene (17ngm~ ). Further, the average concen¬ 
tration of PAHs during different sampling periods 
varied from 15 (Wingene) to 135 ngm“ (Petro¬ 
leumkaai) in the vapour fraction of samples, while 
in the particulate fraction the concentration varied 
from 1 (Mechelen) to 14ngm~ (Zelzate). 

This high concentration of PAHs at Petroleum¬ 
kaai seems to be related with the activities of the 
petroleum industry in the vicinity of this sampling 
site. The sampling site in Borgerhout is located in 
downtown Antwerp nearby a main road with heavy 
traffic. The daily average concentration of PAHs 
reported at this monitoring site shows a similar 
trend for the concentration of individual PAH 
compounds, which observation points to a common 
source of PAHs, most likely the vehicular emission. 
The sampling station at Zelzate (a small town, 
50 km west of Antwerp and the same distance east 
of the North Sea) is primarily considered as an 
industrial site, nearby a main road with a moderate 
traffic compared to the one at the Borgerhout site. 
Samples collected on this site showed great varia¬ 
tions for different sampling days. This may be an 
indication of the influence of fluctuating industrial, 
and/or vehicular emissions combined with the 
variation in the weather conditions. The steel 
industry at Zelzate is located around 700 m south 
of the sampling site and hence the contribution of its 
emission to the total PAH concentration can be very 
fluctuating, mainly depending on the WD and its 
speed. The site selected in Hasselt is relatively close 
to the city centre, and is located at the shore of the 
Albert Channel. The moderately high concentration 
of PAHs at this site may be associated with the 


freighter traffic and their movement at sluice gates 
located next to the sampling site. The sampling site 
in Mechelen (a historic town around 25 km south of 
Antwerp) is of suburban nature and is situated some 
1.5 km of the city centre. The relatively high total 
concentration of PAHs (about the same level as in 
Hasselt) at this site corresponds to the combined 
impact of the industrial influence and the medium 
traffic density. The total average concentration of 
PAHs reported at Wingene is comparatively lower 
than those found at the other sites. This observation 
can be justified, as the sampling site is primarily a 
rural one with negligibly small influence of indus¬ 
trial activities and/or emissions by vehicles, but 
some nearby agricultural activity can contribute to 
the increased background levels of PAHs. 

The average total concentrations of PAHs re- 
ported in the particulate fractions were 5.5 ngnT" , 
while in the vapour fractions it was 48 ngnT" . 
Consequently, most of the PAHs appeared to be 
present in the gas phase (more than 80% of the total 
amount of PAHs). FluT (12%), Chr (11%), BbF 
(11%), BghiP (10%), Ind (8%) and BaP (7%) were 
found to be the most common PAH compounds 
with relatively high concentrations in the particulate 
fraction in Flanders, whereas PA (48%), Flu (14%) 
and FluT (10%) form the major part of the vapour 
fraction. Further, BbF, Chr, BaA and BkF showed 
similar percentages in all filter samples, and hence, 
they seem to originate from a common (type of) 
source at all sites. The PAH compounds with a 
lower molecular weight (i.e., <200 Da), such as PA, 
were found to be almost exclusively in the gas phase. 

3.3. Potential toxic fraction of PAHs in ambient air 
of Flanders 

The carcinogenicity classifications verified by 
EPA Carcinogenicity Risk Assessment Endeavor 
Work Group (EPA, 1994) shows that BaA, BbF, 
BkF, BaP, Chr, DahA and Ind are considered to be 
probable human carcinogens, whereas other PAHs 
such as AcPy, Ant, BghiP, FluT, fluorine, PA and 
Pyr are not classified as promoters of the same 
health risk. Based on these conclusions and our 
observation, it was found that 55% of total PAH 
fractions in aerosols are probable carcinogenic, 
whereas this fraction forms only 2% of total PAHs 
in the vapour phase (Table 3). Hence, it can be 
concluded that most of the probable human 
carcinogens are found to be associated with SPM 
in Flanders. Lighty et al. (2000) have also made 


Table 2 

_o 

Concentration of individual PAH compounds (ngm“ ) at the sampling sites in Flanders during various seasons 
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Table 3 

Potential toxic fraction of total PAH concentrations at various sites in Flanders 


Sampling site 

Phase 

Average PAH concentration 
(ngm -3 ) 

Toxic fraction 3 (%) 

Petroleumkaai 

Aerosol 

3.9 

57.5 


Vapour 

110 

1.5 

Borgerhout 

Aerosol 

8.7 

54.1 


Vapour 

45.7 

1.8 

Zelzate 

Aerosol 

7.4 

61 


Vapour 

44.4 

2.7 

Hasselt 

Aerosol 

4.5 

54.1 


Vapour 

38.2 

1.6 

Wingene 

Aerosol 

0.9 

39.3 


Vapour 

15.6 

2.1 

Mechelen 

Aerosol 

7.4 

49.4 


Vapour 

34.1 

4.3 


a Percent of total PAH concentration. 


similar observation for combustion aerosols. 
Although only 2% fraction of vapour phase PAHs 
is a probable human carcinogen, the health risk 
associated with this portion cannot be ignored as 
the vapour phase PAHs have 10 times higher 
concentration in ambient air in comparison to the 
levels of PAHs associated with aerosol. 

3.4. Influence of meteorological conditions on PAH 
levels 

In order to reveal the atmospheric processes 
responsible for the spatial and temporal distribution 
of air pollutants, it is mandatory to analyse the local 
and the regional meteorology; especially WD and its 
speed, turbulence and atmospheric stability. The 
likelihood of chemical reactions occurring in the 
atmosphere also depends on ambient weather 
conditions because they are influenced by short 
wave radiation, temperature and humidity of air. 
Along with chemical reactions, dispersion and 
dilution processes induce the concentration of 
different substances to vary with regard to time 
and space. Although the relation between meteor¬ 
ology and air pollution has been investigated for 
various pollutants, there is still very little informa¬ 
tion in the literature on the dependence of PAHs 
and urban aerosols (SPM and PM 2 . 5 ) on the 
meteorological conditions and site-specific charac¬ 
teristics. Thus, an attempt is made in this section to 
establish a relationship between these parameters 


through a statistical treatment of the data received 
during the sampling campaigns. 

Table 4 shows the average meteorological condi¬ 
tions during each campaign. The prevailing WDs, 
which may transport moisture or aerosol particles 
from distant sources, play a major role on the daily/ 
seasonal variation of the diverse pollutants. The 
levels of PAHs studied in Flanders also show 
significant variation with relation to WD, although 
analysis of surface WD alone does not adequately 
explain the variability in the concentration of air 
pollutants. Hence correlation analysis was applied 
between various other meteorological parameters. 
This analysis shows that most of the parameters 
bear statistically a significant correlation with each 
other indicating that close association does exist 
between them. Table 5 provides the correlation 
matrix for these meteorological parameters and air 
pollutant levels monitored at different sampling 
sites of Flanders. Based on Table 5, the following 
conclusions can be drawn. 

(I) Particulate PAH levels, as well as the SPM and 
PM2.5, show a negative correlation with tempera¬ 
ture. These features can be explained on the basis of 
mixing height, which reduces with the fall in 
temperature and restricts the mixing of pollutants 
in the atmosphere. Further, more emission can be 
expected during winter from various sources; (II) 
Aerosol and vapour phase PAHs show a negative 
correlation with wind speed, so do SPM and PM2.5. 
Higher WS can lead to less inversion situation and 





Table 4 

Meteorological parameters recorded during each campaign 
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better turbulent dilution of pollutants; (III) The 
particulate PAHs levels develop a highly significant 
correlation with the total SPM and PM 2 5 mass, 
whereas the concentration of PAHs in the vapour 
shows only some correlation with SPM; (IV) The 
concentration of PAHs in vapour and aerosol 
phases also shows some degree of relation and 
likely indicates towards particle/vapour partition¬ 
ing, which mainly depends on the vapour pressure 
of PAHs, amount of fine particles available, 
ambient temperature, etc. 

The correlation matrix for various sites was also 
studied individually, which demonstrates that the 
concentrations of PAHs in the vapour and aerosol 
phases develop a significant correlation at Borgerh- 
out (r = 0.49,/> <0.05), Petroleumkaai (r = —0.56, 
p <0.05), Hasselt (r = 0.68,;? <0.01), Zelzate (r = 
0.92,7? <0.01) and Mechelen (r = 0.41,/><0.05), 
except for Wingene. This is a rural site and probably 
the absence of any significant PAH emission sources 
in and nearby Wingene results in an insignificant 
correlation between vapour and aerosol phase 
PAH levels. At Borgenhout (r = 0.50 ,p <0.05), 
and Petroleumkaai (r = 0.58,7?<0.05), a significant 
correlation was noticed between precipitation and 
RH, probably, because only at these sites some 
showers occurred during the sampling dates and 
periods. At Hasselt, a highly significant negative 
correlation was noticed between WS and the 
concentration of PAHs in both the vapour and 
aerosol phases. Further, at Wingene, the RH shows 
a negative correlation (r = —0.45,;? <0.05) with the 
levels of PAHs in aerosol, while, at Mechelen, the 
RH and precipitation show significant positive 
correlations with the levels of PAHs in the vapour 
phase. Based on these observations, it can be 
suggested that both the site-specific characteristics 
and meteorological parameters play an important 
role in the change of the PAH levels in Flanders. 

3.5. Seasonal variation in the concentration of PAHs 

The concentrations of PAHs during different 
seasons at various sampling sites in Flanders are 
also shown in Table 2. The total PAH concentration 
measured at Petroleumkaai was almost double 
during summer in contrast to winter samples, 
although the concentration of PAHs associated 
with SPM was high in winter compared to the 
vapour phase fraction. The average total concentra¬ 
tions of PAHs at Borgerhout were found to be 
comparable at the beginning and the end of winter. 
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Table 5 

Correlation matrix of various air pollutants and meteorological parameters 



RH 

PPT 

T 

AP 

WS 

A-PAH 

V-PAH 

pm 25 

SPM 

RH 

1.00 

0.15 

-0.10 

-0.04 

-0.29** 

0.03 

0.02 

-0.14 

-0.11 

PPT 


1.00 

0.03 

-0.22* 

0.08 

-0.10 

0.09 

0.10 

0.04 

T 



1.00 

-0.11 

-0.13 

—0.66** 

-0.16 

-0.41* 

-0.43** 

AP 




1.00 

-0.21* 

0.14 

-0.19 

0.17 

0.18 

WS 





1.00 

-0.13 

-0.10 

-0.23* 

-0.23* 

A-PAH 






1.00 

0.20* 

0.62** 

0.63** 

V-PAH 







1.00 

0.19 

0.23* 

PM 2 .5 








1.00 

0.98** 

SPM 









1.00 


**(/?<0.01); *(/?<0.05); n = 92. 

Abbreviations: A-PAH—aerosol phase PAH level; V-PAH—vapour phase PAH level; RH—relative humidity; PPT—precipitation rate; 
T—temperature; AP—air pressure; WS—windspeed. 


There was also no significant change noticed in the 
total concentration of PAHs during different 
seasons at Hasselt and Wingene. 

Other studies conducted during different seasons 
in Europe often showed a higher concentration of 
PAHs in winter than in summer by a factor of 2-5 
(Broddin et ah, 1980; Grimmer et al., 1980; Colmsjo 
et al., 1986; De Raat and De Meijere, 1991). 
Moreover, the concentrations at urban sites were 
found to be 2-3 times greater than at rural locations. 
The present study shows that the concentration of 
PAHs associated with SPM is in good agreement 
with other studies, i.e., showing a higher value in 
winter, while this is not the case for PAHs found in 
the vapour fraction except for Mechelen and Zelzate. 
A possible reason for the high concentration of 
particulate PAHs during winter season seems to be 
related with the very low temperature during these 
periods, which made favourable conditions for the 
condensation/adsorption of these species on sus¬ 
pended particles present in air. 

The vapour phase PAHs seem to be more local in 
origin and hence their levels in ambient air show no 
significant variation during different seasons, 
whereas the PAHs associated with SPM may be 
regional and can be transported far from the source 
of origin. Although some of the meteorological 
parameters (PPT and WS) play an important role in 
controlling the concentrations of PAHs at each site; 
at the same time, the site-specific emission sources 
may also influence their levels in the ambient air. 

3.6. Source identification 

The PAH contents of urban air samples can be 
considered to be a composite of the emission 


products of various local sources combined with a 
contribution from outside areas. Diagnostic ratios 
and PCA are the most widely used techniques to 
identify and characterize the emission sources. 
These mathematical approaches were also proven 
to be important tools in the present study, as 
follows below. 

3.6.1. Evaluation of diagnostic ratios 

The concentrations of specific PAH compounds, 
or a group of PAHs, have been used to indicate 
towards the corresponding emission sources, such 
as Ant, PA, FluT, Pyr, BaA and Chr for coal 
combustion; Ant, PA, BaP and BghiP for coke 
production; PA, FluT and, especially, Pyr for 
incinerators. Ant, PA, FluT and Pyr are indicators 
for combustion of wood; FluT and Pyr for oil 
burning; FluT, Pyr and, especially, BghiP and 
coronene for petrol powered vehicles; FluT, Pyr 
with higher ratio of BbF and BkF for diesel-fuelled 
vehicles (Duval and Friedlander, 1981; Khalili et al., 
1995). These markers indicate some degree of 
similarity and overlap between the profiles from 
different emission source categories. 

Table 6 provides the seven diagnostic ratios 
between total concentrations of PAHs at various 
sampling sites in Flanders. A value of 0.35-0.7 of 
Ind/(Ind + BghiP) ratio has been used for diesel 
emissions (Kavouras et al., 2001). A corresponding 
value has been noticed at all sites with a maximum 
ratio at Hasselt and Wingene. The higher ratio 
(>0.5) of Flu/(Flu + Pyr) has been used for diesel 
emission, while a lower ratio (<0.5) for gasoline 
emission (Rogge et al., 1993; Mandalakis et al., 
2002; Fang et al., 2004). This ratio was found to be 
lower than 0.5 at all sites for vapour phase PAHs, 
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Table 6 

Diagnostic ratios of PAHs at various sites in Flanders 


Location/ 

campaign 

Type of 
sample 

Ind/ 

(Ind + BghiP) 

Flu/(Flu 
+ Pyr) 

BaP / 

(BaP + Chr) 

PA/(PA 
+ Ant) 

BbF/BkF 

BaP/BghiP 

Ind/BghiP 

Petroleumkaai/ 

Filter 

0.49 

0.18 

0.26 

0.92 

2.16 

0.35 

0.94 

1st 


PUF 

0.40 

0.77 

0.48 

0.95 

5.44 

0.07 

0.66 


Total 

0.44 

0.77 

0.29 

0.95 

2.72 

0.18 

0.77 

Petroleumkaai/ 

Filter 

0.33 

0.09 

0.72 

0.95 

2.46 

0.60 

0.48 

2nd 


PUF 

0.44 

0.83 

0.00 

0.95 

2.80 

0.12 

0.79 


Total 

0.34 

0.81 

0.09 

0.95 

2.47 

0.54 

0.52 

Borgerhout/1st 

Filter 

0.30 

0.09 

0.46 

0.78 

1.83 

1.05 

0.42 


PUF 

0.38 

0.66 

0.10 

0.90 

1.71 

0.09 

0.62 


Total 

0.31 

0.65 

0.44 

0.90 

1.83 

0.92 

0.45 

Borgerhout/ 

Filter 

0.31 

0.05 

0.41 

0.89 

2.22 

0.84 

0.46 

2nd 


PUF 

0.48 

0.68 

0.00 

0.91 

0.00 

0.08 

0.92 


Total 

0.32 

0.61 

0.29 

0.91 

2.23 

0.82 

0.46 

Zelzate/lst 

Filter 

0.56 

0.35 

0.40 

0.72 

1.91 

0.91 

1.27 


PUF 

0.45 

0.82 

0.31 

0.93 

2.73 

0.53 

0.83 


Total 

0.54 

0.79 

0.39 

0.93 

2.00 

0.83 

1.17 

Zelzate/2nd 

Filter 

0.41 

0.20 

0.09 

1.00 

2.83 

0.10 

0.69 


PUF 

0.36 

0.66 

0.08 

0.97 

7.00 

0.07 

0.55 


Total 

0.38 

0.65 

0.08 

0.97 

2.97 

0.08 

0.61 

Hasselt/lst 

Filter 

0.65 

0.10 

0.37 

0.91 

2.29 

0.50 

1.84 


PUF 

0.54 

0.79 

0.22 

0.97 

1.29 

0.06 

1.17 


Total 

0.61 

0.77 

0.35 

0.97 

2.21 

0.31 

1.56 

Hasselt/2nd 

Filter 

0.53 

0.02 

0.14 

0.99 

2.48 

0.13 

1.14 


PUF 

0.41 

0.64 

0.00 

0.98 

0.00 

0.00 

0.70 


Total 

0.53 

0.59 

0.11 

0.98 

2.46 

0.13 

1.12 

Wingene/lst 

Filter 

0.61 

0.18 

0.37 

0.84 

2.82 

0.89 

1.55 


PUF 

0.60 

0.80 

0.31 

0.98 

1.83 

0.14 

1.47 


Total 

0.60 

0.78 

0.35 

0.98 

2.65 

0.40 

1.50 

Wingene/2nd 

Filter 

0.42 

0.54 

0.11 

0.98 

2.75 

0.06 

0.72 


PUF 

0.32 

0.73 

0.01 

0.98 

— 

0.07 

0.47 


Total 

0.39 

0.71 

0.02 

0.98 

3.00 

0.06 

0.64 

Mechelen/1 st 

Filter 

0.43 

0.45 

0.11 

1.00 

4.04 

0.13 

0.76 


PUF 

0.37 

0.80 

0.03 

1.00 

1.80 

0.20 

0.59 


Total 

0.41 

0.79 

0.04 

1.00 

3.46 

0.16 

0.69 

Mechelen/2nd 

Filter 

0.32 

0.04 

0.41 

0.91 

2.21 

0.72 

0.48 


PUF 

0.43 

0.69 

0.00 

0.92 

4.00 

0.15 

0.75 


Total 

0.33 

0.58 

0.25 

0.92 

2.22 

0.70 

0.49 


likely indicating their origin from gasoline cars, 
while the concentration of PAHs associated with 
particulate matter shows its origin from diesel 
emissions. Khalili et al. (1995) and Guo et al. 
(2003) indicated that the ratio of BaP/(BaP + Chr) 
was 0.49 and 0.73 for diesel and gasoline engines, 
respectively. Our results do not correspond to these 
values and indicate that other pollution sources 
possibly interfere. A higher ratio of BbP/BkP has 
been noticed at all sites, which indicates the presence 
of diesel emission. The BaP/BghiP ratio higher 


than 0.6 refers to the presence of traffic emission 
and contribution from other PAH sources (Pandey 
et al., 1999; Park et al., 2002). A similar value of this 
ratio was found only at Borgerhout, while the other 
sites in Flanders do not match to this data. 
Caricchia et al. (1999) reported that the Ind/BghiP 
ratio for gasoline engines is about 0.4, while the 
ratio for diesel engines approaches 1.0. Most of the 
values for this ratio were higher than 0.5, and this is 
probably in agreement with the diesel emission. 
Although it is clear from the study of diagnostic 
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ratios that the diesel and gasoline emissions are 
significant emission sources of PAHs in Flanders, it 
also indicates that other pollution sources exist and 
contribute to the PAH levels as well. Hence, in 
order to enhance the accuracy of the emission 
source identification, the method of PCA was 
applied as follows in the next section. 

3.6.2. Principal component analysis 

Using PCA, it is possible to simplify the inter¬ 
pretation of complex systems and to reduce the set 
of variables to few new ones, called factors. Each of 
these factors can be identified as either an emission 
source, or a chemical interaction. Many of these 
factors indicate more than one possible cause. In 
this study, each factor from PCA is associated with 
a source characterization by its most representative 
PAH compound(s). Source groupings were deter¬ 
mined by using PCA with varimax rotation and the 
retention of principal components having eigenva¬ 
lues higher than 1 of the complete data set of PAH 
concentrations. Further, the principal components 
that showed more than 15% of total variance of the 
data set were used as factors. Loadings determined 
the most representative PAH compound in each 
factor and a value of higher than 0.5 was selected. 

Table 7 provides the results of factor analysis of 
total concentrations of PAHs at all the sites studied 
in Flanders. The PCA results show that two of the 
factors explain the main part of the data variance; 
therefore, one element for each factor has been 
chosen as a tracer. The high factor loading of FluT, 
Pyr and especially of BghiP and coronene has been 
suggested for gasoline powered vehicles (Duval and 
Friedlander, 1981; Masclet et al., 1986; Khalili 
et al., 1995), whereas Guo et al. (2003) also include 
BaA, BaP, BbF, BghiP and Ind as source markers 
of gasoline emission. Diesel emission has high factor 
loading for FluT, PA, Ant and Pyr (Caricchia et al., 
1999; Yang et al., 1998; Ho et al., 2002; Omar et al., 
2002; Fang et al., 2004). Other studies also suggest 
that FluT, Pyr with higher factor loading of BbF 
and BkF indicated towards diesel-powered vehicles 
(Duval and Friedlander, 1981; Khalili et al., 1995). 
Based on these observations, the factor 1 of filter 
samples at Borgerhout, Hasselt, Mechelen and 
Wingene indicates towards vehicular emission. 
Further, the factor loadings of Petroleumkaai and 
Zelzate also seem partly related with vehicular 
emission. 

According to Daisey et al. (1986), volatile 
compounds such as FluT and species of high 


molecular weight such as Ind, probably occur 
together in the combustion of lubricating oil, which 
might be emitted by industrial sources. The high 
factor loading of AcPy, Acp, Chr, BbF and 
benzo(e)pyrene, indicates towards stationary emis¬ 
sion sources (Yang et al., 1998; Kulkarni and 
Venkataraman, 2000). The factor 2 of filter samples 
seems to show a high factor loading for wood/coal 
emission at Borgerhout and Hasselt, whereas in 
Mechelen it indicates towards incinerator/oil burn¬ 
ing, in Petroleumkaai towards oil burning and in 
Wingene towards stationary emission sources. 

The high level of Ant, PA, BaP, BghiP and Chr 
has been suggested for steel industry emissions. The 
PAHs compounds from coke ovens are released 
during coking operation as fugitive emission in steel 
industries. The observation of high factor loading of 
selected PAHs at Zelzate explains that the levels of 
PAHs have a relation with local activities in 
surrounding areas, i.e., the presence of steel industry 
and vehicular traffic. Further, the high factor 
loading of related compounds in vapour phase 
samples also confirm this assumption, whereas the 
factor 2 of filter and PUF samples also indicates 
towards secondary emission sources, i.e., oil burn¬ 
ing and stationary sources. 

The Factor 1 for PUF samples indicated towards 
diesel (Borgerhout, Mechelen), gasoline (Hasselt), 
petroleum and oil burning (Petroleumkaai) and 
vehicular emissions (Wingene and Zelzate) as 
primary sources in Flanders. Whilst factor 2 for 
PUF samples suggests gasoline (Borgenhout, Me¬ 
chelen), wood/coal (Hasselt), stationary (Wingene, 
Zelzate) also occur as secondary sources for PAHs 
emission. 

In general, the source marker study for both 
vapour and particulate phase PAHs reveals that 
vehicular emission seems to be dominating over 
other sources of PAHs in Flanders. In most of the 
cases, vapour phase PAHs represent more local 
source of origin and can also further specify the 
particular dominating sources, i.e., gasoline or 
diesel. Most of the studies, conducted to identify 
the potential sources of PAHs, have used total PAH 
levels for PCA, which likely give similar results, but 
it is suggested to use the application of PCA 
separately for both vapour and particulate phase 
PAHs to have more clear indication about the local 
and regional emission sources of these compounds. 
This consideration is also true for the statistical 
analysis of other material systems existing in more 
than one phase in the atmosphere. 


Table 7 

Factor analysis of total PAH concentration data set of various sampling sites 
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4. Conclusion 

Samples collected from various sites in Flanders 
were extracted using PLE prior to their analysis 
with HPLC. This method offers a very fast 
extraction of PAHs from the air samples with 
excellent recovery and minimum solvent consump¬ 
tion. Fast procedure for PAH extraction and least 
consumption of solvent with equal or more extrac¬ 
tion efficiency to conventional methods offer PLE 
as a fruitful choice of extraction. 

The concentration of PAHs at various sites in 
Flanders was found to be related with site-specific 
characteristics. In general, both the Petroleumkaai 
and the Borgerhout sampling sites showed a 
relatively high level of PAHs. Additionally, the 
concentrations of PAHs were found to be signifi¬ 
cantly higher (approximately 10 -fold) in the vapour 
fraction of samples than in the particulate fraction. 
This means a certain raise in the risk of potential 
exposure of human to these compounds, in addition 
to the health effects of PAH fractions associated 
with SPM and PM 2 . 5 . It is to be noted here that the 
vapour phase PAHs well describe the site-specific 
activities and hence they can be used as more 
reliable indicators of the local PAH emissions. 
FluT, Chr, BbF, BghiP, Ind and BaP were found 
to be the most common PAH compounds with 
relatively high concentrations in the particulate 
fraction, whereas PA and FluT form the major 
part of the vapour fraction. The emission from 
vehicle exhausts (diesel, gasoline) was found to be a 
significant contributor to the atmospheric levels of 
PAHs in the urban environment. 
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